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We propose a mean-field method to calculate approximately the spacing distribution functions 
p( n \s) in ID classical many-particle systems. We compare our method with two other commonly 
used methods, the independent interval approximation (II A) and the extended Wigner surmise 
(EWS). In our mean- field approach, p^ n \s) is calculated from a set Langevin equations which are 
decoupled by using a mean- field approximation. We found that in spite of its simplicity, the mean- 
field approximation provides good results in several systems. We offer many examples in which the 
three methods mentioned previously give a reasonable description of the statistical behavior of the 
system. The physical interpretation of each method is also discussed. 
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I. INTRODUCTION 

The spacing distribution functions p( n \s) are often 
used to describe the statistical behavior of many-particle 
systems in one dimension (ID) [1-7J. By definition, 
p( n \S) is the probability density that an interval of 
length S which starts at a particle contains exactly n 
particles and that the next, the (n + l) th particle, is in 
[S, S + dS]. The relative spacing is defined as s = S/ (S), 
where (S) is the average of S. However, p( n \s) are 
also used in other contexts. For example, in ID sys- 
tems with domains, S represents the spacing between 
domain boundaries [SHE], while in random matrix the- 
ory and in quantum systems this variable represents the 
spacing between adjacent energy eigenvalues [16-21]. In 
the physics of surfaces, S can be the distance between 
islands in epitaxial growth models or the terrace width 
between adjacent steps on vicinal (misoriented) surfaces 

I22H25]. 

The spacing distribution functions are useful even in 
non-equilibrium systems having dynamical scaling. One 
such system is the coalescing random walk (CRW). In 
CRW all particles execute independent random walks, 
suffering a fusion reaction (A + A — >• A) when two parti- 
cles meet. Clearly in this system the number of particles 
and (S(t)) are time dependent. Despite this, it is possible 
to define a scaled spacing distribution according to 



pM(8) = (S)pW(8(S),t). 



(1) 



Since p^ (s) in Eq. (yj) does not depend on t, it can be 
compared with the spacing distribution functions of an 
equilibrium system. More information about the CRW 
is given in Ref. [26] . 

In this paper we propose a mean-field model to calcu- 
late in an approximate way the spacing distribution func- 
tions of classical many-particle systems. Our method is 



compared with other existing methods commonly used 
to calculate p( n \s). These approximate methods are im- 
portant because they are often used to obtain the spacing 
distribution functions for systems where the exact analyt- 
ical solution either cannot be obtained or, if it is obtained, 
cannot be handled easily because of its complexity. In 
Section II, we provide a brief introduction to the spacing 
distribution functions, giving some important definitions. 
In Sections III and IV we briefly review the independent 
interval approximation (IIA) and the extended Wigner 
surmise (EWS) [27 , respectively. (NB, the EWS dif- 
fers [27] from the so-called generalized Wigner surmise 
(GWS) [22-24 , which is not used in this paper.) In Sec- 
tion V we develop two mean- field models. In each case 
we provide several examples of the use of these methods, 
analyzing their advantages and limitations. Finally, in 
Section VI we offer conclusions. 

All our numerical data for Dyson's Brownian motion 
model were generated with a Monte Carlo simulation that 
uses the standard Metropolis algorithm. We used a lat- 
tice with L = 10000 sites with N = 50 particles. The 
statistics take into account over 20000 realizations. 



II. BASIC DEFINITIONS 

Consider N particles which can move around a cir- 
cle of circumference L; periodic boundary conditions are 
imposed, that is, xn+j = Xj, where xj is the position 
of the j-th particle. If the system is in equilibrium at 
an inverse temperature /3, then its statistical behavior 
is totally defined by the joint probability distribution 
P N (xi, • • • ,x N \fS) 



P N (xi, • • • ,x N \0) 



1 



Z N (L; i 



-(3 VOi,--- ,x N ) 



(2) 
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where V(xi, — ' , x/v) is the total interaction energy 
among the N particles, and Zn is the conflgurational par- 
tition function of the system. Then Pjy (#!,••• ,xjsr]P) 
represents the probability density to find particle 1 in 



2 



[xi, x\ + dxi], particle 2 in [x2, %2 + dx^ etc. In prac- 
tice this joint probability distribution cannot be obtained 
easily from experiments or numerical simulations because 
it depends on many variables. Instead, one usually eval- 
uates the spacing distribution functions p( n \S), where 
n > 0. Each p( n \S) contains reduced information about 
the system. In order to have a more complete description 
of the statistical behavior of the system, it is necessary 
to know p( n \S) for all n. From them it is possible to 
calculate the pair correlation function g(s) by using the 
expression 



300 



p( n )(s). 



(3) 



n=0 



The sum in Eq. ^ seems to be a formidable task. How- 
ever, in many cases g(s) quickly relaxes to 1 and there- 
fore, just the first p( n \s) have to be calculated explic- 
itly. Many thermodynamic quantities can be expressed 
in terms of the pair correlation function and the interac- 
tion potential between particles [28] [29] . 

When the N particles interact via a pair potential v(r), 
the total energy of interaction reduces to 



N q 



V(x lr - ,x N ) = 2^ 2_^v(x m+j 



(4) 



with q the number of interacting neighbors. For q = 1 
we have nearest-neighbor interactions and for q = TV — 1 
each particle interacts with all other particles. We hence- 
forth denote the latter as "full-range" interactions, cor- 
responding to infinite range in the thermodynamic limit. 
As mentioned previously, the probability density to find 
the N particles around the positions x\, • • -xn is given 
by Eq. As usual, Zn can be calculated from 

Zn(L; (3) = dx\ ■ ■ ■ dxjy S (A) P/v (#i, • • • , xn] /?) , 



(5) 

where A = L — ^2 i=1 — Xi)- Making the change of 

variables Si = Xi+i — xi, the partition function takes the 
form 

Z N (L; /3) = y d5i • • • dS N 6 (A) P N (Si, • • •, S N ; (3) , 



( 6 ) 

now with A = L — ^2 i=1 Si. In the same way, the joint 
probability distribution can be written as 



with 



N 



= H S m)+v(S m + S m+ l) 



m=l 



■v(S n 



■Sm+ q -i)]. (8) 



Taking f(S]/3) = f(S) = e~^< s \ we find 

N 



Pn (Si, • • • , Sn', f$) 



1 



Y[ F(Srn, "' , S, 



m+q-l) 



(9) 

where we define F(S m , • • • , S m + q -i) f(Sm)f(S m + 

Sm+i) ■ • • f(S m H 1- S m + q -i). The joint probability 

distribution of n consecutive spacings P n (Si, • • • , S n ; /3) 
is given by 

Pn (Si, • • • , S n ; f3) = J dS n +i • • • dSjsj Pn (Si, • • • , Sat; 0) . 

(10) 

By definition, the n th spacing distribution function 
p( n \S) = p^ q \S) can be written as 

poo 

p(»i«)(5)= / dS 1 ---dS n+1 S( V )P n+1 (S 1 ,--- ,S n+1 ;f3), 
Jo 

(11) 

with r] = S — Si- Note that this notation makes ex- 

plicit the dependence of the spacing distributions on the 
number of interacting neighbors q. The average spac- 
ing between particles is (S) = L/N. Then, the scaled 
probability density is 



^ q \s) = 

Jo 



dSi---dS n+ iS(X)P n+ i (Si,--- i S n +i \ / 



(12) 

with A = r)l (S). Note that Eq. (12) satisfies the normal- 
ization conditions [17] 



dsp^ q) (s) = 1 and / dssp^ q) (s) = n + l. 
o Jo 

(13) 



In general, the integral given in Eq. (12) cannot be 



easily calculated by analytical methods. In fact, it can be 
solved just in few cases. A simple example of this is given 
by the potential v(r) = — ln(r) with q = 1. By using the 
Laplace transform method it can be shown [TJ [18] that 



\(l+0)(l+n) 



,£+n(l+/3) -s(l+/3) 



r[(l + /3)(l + n)r ^ " (M) 

The case q — N — 1 — > oo was solved by Dyson and so is 
termed "Dyson's Brownian Motion Model" [30], [31] . The 
general case, 1 < q < N — 1, with AT, L — >• oo was solved 
in Ref. [1 through integral equations. Henceforth, we de- 
note as Dyson's Brownian model a ID system in which 
the particles interact through a logarithmic potential, re- 
gardless of the range of interaction q. 



III. INDEPENDENT INTERVAL 
APPROXIMATION 

A. General Formalism 

The easiest method to obtain approximate expres- 
sions for p( n]q \s) is the independent interval approxima- 
tion (HA). In the IIA, all spacing distributions are gen- 
erated from the nearest-neighbor distribution p^ q \s). 
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The essence of the approach is to neglect the correla- 
tions between the sizes of adjacent spacings. Conse- 
quently, the probability to find particles around positions 
x\ , X2 , • • • 5 x _/v is given by 



N 



Pf(x u x 2r .- ,x N ;/3) = Y[p^ q) (x i+1 - Xi ) 



2=1 



~P Z)iLl V e f{(x i+1 -Xi) 



(15) 



Comparing Eq. ( |15| ) with the Boltzmann factor given in 
Eq. ([2|, we conclude that 



N 

V(xi, • • • ,X N ) + C = ^V e ff(^i+l - Xi). (16) 

2=1 

The constant C is needed to ensure normalization of the 
joint probability distribution. 

In IIA, each particle is taken to interact with its near- 
est neighbors via an effective pair potential given by 
^eff(^) = — In (p( 0;gr )(r)). This effective potential is gen- 
erally different from the real potential; most importantly, 
the effective potential usually depends on the inverse 
temperature of the studied system. Note that the IIA 
reproduces exactly the functional form of the interac- 
tion potential V(#i, ■ • ■ ,%n) m the case of q = 1. How- 
ever, we caution that entropic repulsions (arising from 
fermionic non-crossing interactions su ch as occur in the 
TSK model of steps discussed in Sec. VC) intrinsically 
involve all particles; hence, this part of V(x±, • • • , xn) in 
Eq. (16) cannot be decomposed exactly into components 



between neighboring parts. 



B. Applications 

As a sample application of the IIA, consider Dyson's 
Brownian model with q = 1. From Eq. ( [l2| ) the statistical 
behavior of this model is clearly equivalent to that in 
which particles interact only with their nearest neighbors 
through an effective pair potential given by 



v e s(r) = -ln(r) + 



C, 



(17) 



In Eq. (17), the first term represents the repulsive force 



between adjacent particles while the second takes into 
account the effect of the average pressure force on the 
ends of the interval of length r due to the particles which 
are outside of the interval. The statistical behavior of 
Dyson's Brownian model with nearest-neighbor interac- 
tions in the IIA is given by 



N 



Pf(x u x 2 , 



,x N ;f3)=C N Y[< 



,/31n(S J -)-(l+/3)S J - 



(18) 



From this discussion it should be clear that many systems 
can share the same statistical behavior even when they 



have different interaction potentials between particles. In 
the above example, the real and effective system have the 
same statistical behavior but their interaction potentials 
differ. 

As for all systems with just nearest-neighbor interac- 
tions, in the IIA all spacing distribution functions can be 
calculated analytically in Laplace space [TJ [2] . This ap- 
proximation should work properly in systems with short- 
range interactions between particles because the correla- 
tions between the sizes of adjacent spacings are usually 
weak compared to the case of long- or infinite-range in- 
teractions. 

Consider again Dyson's Brownian motion model. By 
integrating Eq. (18) we see that, for q = 1, the IIA de- 



scribes exactly the statistical behavior of this system. 
Fig. [TJa) shows that the IIA is also a good approxima- 
tion for q = 2. The functions plotted in Fig. Ilia) were 



calculated from the exact expression for p^ 0;2 \s) given m 
Ref. pQ for q = 2 



p (0;2) (s) = s (2.4773 + 6.0681 s + 3.7159 s 2 ) e 



-3s 



(19) 



through Eqs. ([12| and (15). In this case, the exact ex- 
pression for p^'^s) is given by 



p (1;2) (s) = s A (2.5054 + 3.068 s + 0.7516 s 2 



while in the IIA we have as result 



-3s 



> (20) 



,(i;2) 



s 3 (1.02247 + 2.50462s + 2.14728s 2 + 
0.751391s 3 + 0.0986s 4 ) e~ 3s . (21) 



3s 



We conclude that the IIA reproduces the p( 1;2 ^ (s) ex e 
behavior of the spacing distribution functions for large 
values of s but fails for small values of s, predicting 
p( 1;2 )(s) oc s 3 rather than the exact result p( 1;2 )(s) oc s 4 . 

However, as q increases the IIA becomes very poor, as 
shown in Fig.[TJb) where we took the widely used Wigner 
surmise expression for q = N — 1 — >> oo 



p(0;oo) (s) = | se 



(22) 



and calculating p( n '°°\s) from Eqs. (12) and (15). As 



expected from the lack of correlation, the IIA becomes a 
progressively poor approximation as n increases. 

Another example of a system in which the IIA can be 
applied successfully is the point island model for epitaxial 
growth [25] I32U35] . In this model, the monomers are 
deposited at random on a d- dimensional lattice. The 
monomers diffuse across the lattice until reaching a site 
occupied by another monomer, at which point an island 
is formed. This coalescing process is called nucleation. 
While the monomers are mobile, the islands remain static 
and do not grow laterally. Fig. [2] shows that the IIA is 
a good approximation for this model. In Fig. |2|a) we 
use the heuristic equation p^°°\s) = A a s a e~ BaS with 
a = 3/2 given in Ref. [33] while in Fig. ^b) we use 
the same expression but with a = 2 [36 j. Additional 
examples of the IIA can be found in Refs. [I2j [TSJ [37] . 
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FIG. 1: (Color online) The II A for Dyson's Brownian motion 
model with a) q = 2 and b) q = 49. 



In short, in the II A the statistical behavior of an 
iV-particle system is modeled by using another system 
whose particles interact through an effective potential 
with their nearest neighbors. The effective potential is 
a function of the inverse temperature and must be cal- 
culated for each particular system. However, as shown 
in Refs. [12] [15], p(°' ,q \s) does not uniquely define the 
system. 



IV. THE WIGNER SURMISE 

The exact expression for the spacing distribution func- 
tions for Dyson's Brownian motion model with full- 
range interactions is complicated and unwieldy. How- 
ever, Wigner proposed accurate and simple expressions 
for p( 0;oo )(s) in the special cases of f3 = 1, 2 and 4. Abul- 
Magd and Simbel [17] subsequently extended Wigner's 
surmise to n > by making the ansatz 



P 



(n+l;oo) 



OC S 



Jo 



9 (n;oo) (s), (23) 



FIG. 2: (Color online) The IIA for the point island model of 
epitaxial growth, d = 1 on the top and d = 2 on the bottom. 
We also show the EWS approximation for this model in the 
2D case. 



using Wigner's expressions for p( 0;oo \s), and assuming 
Gaussian decay. In this EWS, p( n] °°\s) for Dyson's 
Brownian model with q = N — 1 — oo is written as 



P 



(24) 



with a n = n+(n+l)(n+2)/3/2. The constants A n and B n 
are calculated from the normalization conditions given 
by Eq. O. In the EWS, the functions p^°°\s) are 



completely described by their behavior in the two lim- 
its 8^0, where p( n;oo \s) oc s ftn , and s -+ oo, where 

As mentioned above, the EWS was developed to find 
analytical expressions for the spacing distributions of 
Dyson's Brownian motion model. Nonetheless, it can 
also be applied to a great variety of systems with the 
suitable choice of a n . One advantage of the EWS com- 
pared to alternatives like the IIA lies in its mathematical 
simplicity. Another is that the IIA is well defined only 
for ID systems while the EWS does not make any as- 
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sumptions about the geometry of the system and it can 
be applied to more general systems. In Fig.[2|b) we have 
applied the IIA simply as a mathematical artifact for a 
2D case; in this case it is not possible to define spacings 
between islands. 

Fig. [2^b) shows that the Wigner surmise is a good 
approximation for the point island model for epitaxial 
growth with d = 2. Additionally, in many cases the 
EWS captures the essential physics of the system. E.g., 
for Dyson's Brownian motion model with infinite-range 
interactions, the EWS for n = is accurate to better 
than of the 5% for f3 = 1 and even better for j3 — 2 or 4 
[38] . For n > the fits visually seems to have the same 
quality as n = [12] [T7]. In the IIA all p^ q \s) can 
be calculated in Laplace space, but the inverse of those 
expressions usually cannot be calculated for large values 
of n. For the point-island model (d = 2) in the EWS, 
we invoke Eq. (24) as the ansatz and calculate a n from 



a simple argument: Since we know that the IIA gives a 
good fit for this model, we can expect that for s <C 1 



(n;oo) 



poo poo n _^~} 

/ d Sl --- ds n+1 TTp("'°°)( Sl )cx S 3 ™ +2 ) 

Jo Jo i=1 

(25) 

so that a n = 3 n + 2 [36 . In Fig. [2^b) we see that the 
differences between the IIA and the EWS are almost im- 
perceptible for the point island model in d = 2. 

Unfortunately, the EWS is not always a good approx- 
imation. For example, consider the CRW and Dyson's 
Brownian model with infinite range of interactions and 
(3 = 1. In Ref. [39] it was shown that both systems have 
Gaussian behavior for large values of s and that a n is 
given in Eq. (24), with j3 = 1, for both systems. How- 



ever, the EWS gives good results for Dyson's Brownian 
model but not for the CRW. 

Another example is Dyson's Brownian motion with 1 < 
q < N — 1, where p( n]q \s) ex e~ BnS for s — » oo pQ. 
It is then tempting to extend our previous approach by 
writing p^ n ' q \s) in the more general form 



V 



(n;q) 



(s) = A n s a "e- B " s \ 



(26) 



here with 7 = 1, while a n is unchanged from Eq. (24). In 
fact, for finite-range interactions, the functions ^'(s) 
have the form 



p^"\s) = A n s a "Q n (s)e- 



(27) 



where Q n (s) is a polynomial whose degree depends on 
n and /3 [1 . Thus, for Dyson's Brownian model with 
finite-range interactions, the spacing distribution func- 
tions cannot be described solely by their behavior in the 
limits 8—^0 and s —> oo; the behavior of p( n;q \s) for 
intermediate values of s is clearly important. 

Additional applications of the EWS can be found in 
Ref. [39 , where the EWS was applied to study the sta- 
tistical behavior of an N particles system with nearest- 
neighbor interaction v(r) = ^r 2 — In (§r) . In this case 
the EWS gives excellent results with the ansatz 7 = 2 



and a n = 2 n + 1. This potential was used to model two 
d = 1 non-equilibrium systems with domain formation 



V. MEAN-FIELD MODEL 

In Refs. [22] [23] a mean-field (MF) model is proposed 
to obtain an analytical expression for the distribution of 
terrace widths (spacings between adjacent steps) on vici- 
nal surfaces. Following these ideas, we propose a new ap- 
proximate description for Dyson's Brownian model with 
finite-range interactions, that is, for arbitrary values of 
q. This system was solved exactly in Ref. [1 . However, 
the expressions found there for the spacing distribution 
function in the range 1 < q < N are unwieldy because 
they involve integral equations which increase in diffi- 
culty for large values of q and /3. Rather (and following 
Refs. [22j[23]), we propose a L mean-fieW approximation 
to obtain simple expressions for this model. 



A. Nearest-Neighbor Distribution with Finite 
Range Interaction 

1. Simple mean-field approximation (SMF) 

According to the logarithmic potential proposed by 
Dyson, the force F m = F m (xi, • • • ,xn) on the m th par- 
ticle is given by 



E 



1 



1 



(28) 



Consequently, the deterministic damped-oscillator equa- 
tions for the m th particle are 



M 



d 2 a 



2 dX nr 



dt 2 



r dt 



(29) 



where 2 /Y — anticipating the subsequent development — is 
the friction coefficient. In the overdamped limit of strong 
friction, the system relaxes quickly to the equilibrium, 
allowing us to take d 2 x m /dt 2 = and rewrite Eq. (29) 
as first-order differential equations 



dX nr. 



dt 



(30) 



As in Ref. [40] , we then consider the time evolution of the 
spacings between adjacent particles, S m = x m+ i — x m , 
adding a stochastic term to produce the set of Langevin 
equations 

-^ = - { F m+1 -F m ) + - Vm (t), (31) 

where the r) m (t) introduce Gaussian white noise for each 
spacing, and 
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m+l 



ra+i 



^2i=l Sri 



(32) 



Taking the average () in Eq. (31) as discussed in the 
Appendix [Aj we have 



dS 
~dt 



5 + (j-l)<5>' 



(33) 



where the Cj's are constants which can be interpreted as 
the renormalization of the strength of interaction. The 
first sum at the right of Eq. ( 33 ) can be calculated easily 



dS 
~dt ~ 

where H„ 



CiT 

(S) 



H p -c 2 rj2 



iS + v-i)(S)' 



(34) 



_1 is the generalized harmonic 



number of order q. After the change of variables r = 
Yt/ (S) 2 and s = S/ (5), it is straightforward to show 



ds 



(35) 



Then the effective potential associated with Eq. ( 35 ) is 

q 

Ve$(s) = C 1 H q s-C 2 ^ln( S + (j-l)) + C 3 , (36) 
i=i 

The nearest-neighbor distribution ^^(s) is given by 
q 



p^\s) = A l[( s + (j-l)) 



/3C 2 e -pdH q s_ 



(37) 



The constant C2 can be calculated by using the fact 
that for low densities p( 0;q \s) « e _ ^ v ^ s \ which leads to 
(72 = 1. The remaining constants can be calculated by 
using the normalization conditions. In the SMF model, 
Eq. (37) has a similar functional form to Eq. (|27|). How- 



ever, the order of the polynomial factor is not the same 
in both cases. For example, for the case of q = 2 and 
/3 = 1, the SMF model gives 



p^ 2 \s)^ 3.30 s(l + s) e- 2 - 45s , 
which clearly differs from Eq. (Tl9j) . 



(38) 



2. Improved mean- field approximation (IMF) 

We expect that for large values of q the SMF approx- 
imation begins to fail. The reason lies in the fact that 
implicit in Eq. (32) there are more terms which depend 



on S m . For example, in the case q = 2 we have such 
terms as 



1 



1 



%m-\-l 



(39) 

which clearly depends on S m . The number of such terms 
increases with q. Furthermore, for full interactions be- 
tween all pairs of particles, Eq. (37) does not reproduce 
the behavior of the Wigner surmise. Consequently, we 
must formulate an improved mean-field approximation 
(IMF hereafter). After some algebra, it is possible to 
write Eq. (31) as 



dS m 
~dt 



■p 9-1 

2 . =1 {x rn j r jj r i 

n q-1 



«^m+l ) (^m+j + 1 



*^m) 



\ fa, 



^m—j ) («^m+l 



1 



1 



r r 



»?m(*)- 



(40) 



Taking the average () in Eq. (40) as discussed in the 
Appendix [A} we find 



dS 
~dt 



9-1 



r E7 



Ci5 



rci rc 2 



5 ■ 



(41) 



where the constants C\ and C2 are defined implicitly in 
Appendix [aJ After the changes of variables r = Yt/ (S) 2 

and a = S/ (S), then using Y%Z\tiU + 1)] _1 = (<?" !)/<?, 
it is straightforward to show 

(7-1 



(is 



s 



^1 



s 



77. 



(42) 



Then, the effective potential associated with Eq. (42) is 

2 



1>eff(s) 



Ci 



C 2 ln (s) + f (^)f )+ ^s + C. (43) 
it is easy to find C 2 = 1. 



Again, using p( 0;9) (s) « e 

The nearest-neighbor distribution £>( 0;9 )(s) is given by 



p(° ; «)(s) = As^exp 



C x (q-\ 



(s 2 ) 



As 



(44) 
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In the case of full-range interactions, the nearest- 
neighbor distribution can be written as 



P 



(0;<z) 



(s) = As p e 



P-{PC 1 /2){s 2 /{s 2 )) 



(45) 



By using the normalization conditions, one can show that 
Eq. (45) is equivalent to the Wigner surmise described 
by Eq. (24) for n = 0. In the case of nearest-neighbor 

. For 

> 



q 



interactions, Eq. (44) takes the form of Eq. ([14]) . 
1, Eq. (44) predicts a Gaussian tail for 



arge val- 
ues of s instead of the exponential in Eq. ([271). Despite 



this discrepancy with the exact result, it is reasonable to 
expect that the IMF model given by Eq. (44) leads to 
good results, especially for large and small values of q 
(where it reduces to the Wigner surmise and to Eq. (14), 
respectively). 

Before generalizing our mean- field models to n > 0, we 
offer some comments about Eq. (44). First, as mentioned 
before, the IMF for the case of full-range interactions 
was used previously [22 , 23 to describe the terrace-width 
distribution between steps on vicinal surfaces. However, 
Eq. (44) corresponds to a generalization for the case of 



finite-range interactions. 

Second, the capture zone (CZ) distribution of islands 
generated in the early stages of epitaxial growth was 
modeled excellently by using the functional form given 
by Eq. (44) [41 . There, the authors use the maximum 



entropy method to justify this functional form as an ap- 
proximation to the CZ distribution. 

Finally, Eq. (44) has the same functional form as 
Eq. (4.1) in Ref.^42 , where this expression is proposed 
for Dyson's Brownian model with < j3 < 4 based on 
heuristic arguments and on the asymptotic form found by 
Dyson for p^°°\s) [31] • They found that Eq. ^ gives 
good results for Dyson's Brownian model with complete 
range of interactions for the values of j3 mentioned pre- 
viously [42]. From our results, this asymptotic form can 
clearly be interpreted as a mean-field approximation for 
Dyson's Brownian motion for large values of q. 



We emphasize that Eq. (44) is a generalization of the 



Wigner surmise for finite-range interactions, which seems 
to be related to many different models. 



n < q — 1 and n > q — 1. For the first we find 
_ /3(n + l)(n + 2) , 



(47) 



which is the same exponent as for random-matrix ensem- 
bles and Dyson's Brownian model with q — > oo and /3 =1, 
2, or 4. However, for n > q — 1 we find 



f3q(3 + 2n-q) 
a n = h n, (48) 



differing from its counterpart in random matrices. Hence, 
the SMF model for arbitrary values of n can be written 

as 

p(^)(s)=As a "e-e c ? H « s ]J(s + jf (49) 

while the IMF model is given by 

p^( S ) = A n s a "e~^^ ) ^:~^ S (50) 



where a n is defined by Eqs. (47) and (48). The results 
of the SMF and IMF models for j3 = 1, 2 and 4 are 
shown in Figs. [3] [4] and [5] respectively. In general, the 
IMF model gives better results than the SMF model. 
However, for large n and small q both give similar results 
(see Figs. |4^a) and [5|. The SMF model describes the 
spacing distribution functions reasonably well for small 
values of q and j3 (see Figs.[3^a),[3jb) and[4^a)). However, 
as expected from discussion at the beginning in Section 
V A2[ the SMF model gives poor results for higher values 
of q and /?, as shown in Fig. [3jc). In general, the IMF 
gives good results. For large values of g, the IMF model 
gives excellent results, even for large values of n, at least 
for P = 1, 2 and 4. However, for small values of q and 
(3 = 2,4, Figs. |4^a) and Figs. [5|a) show that it gives 
good results only for small n. The SMF and IMF models 
give the exact statistical behavior of the system in the 
case q = 1 while for full-range interactions, only the IMF 
reduces to the EWS. 



B. Arbitrary Spacing Distribution Functions 



For arbitrary interaction range, the analytic expression 
for the n th spacing distribution function is not so clear. 
The spacing distribution functions in the limit s — 
have the form 

p^ q) (s)(xs an , (46) 

where a n is a function of /3, n and q. From Eq. ( fl2| ) it 
is possible to find this exponent by following the same 
kind of arguments used in Ref. [39] for the CRW case, 
see Appendix [B] However, for Dyson's Brownian motion 
model it is necessary to distinguish between the two cases 



C. Another Example: TSK model 

Vicinal crystals usually have terraces oriented in the 
high-symmetry direction separated by steps of atomic 
height, see for example Ref. [43]. The terrace-width dis- 
tribution (TWD) of vicinal surfaces, that is, the distri- 
bution of separations between adjacent steps has special 
interest for experimentalist because it gives information 
about the the interaction between steps [24]. It is as- 
sumed that the steps cannot cross each other and, in the 
simplest case of the terrace-step-kink (TSK) model, that 
the dominant excitation is the formation of kinks. Typ- 
ically the step edges interact via a pair potential v(x), 
where x is the distance between the steps. Usually, the 
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(b) 



(c) 




(d) 




FIG. 3: (Color online) Comparison of SMF and IMF models with simulated data, for f3=l, for interactions spanning the range 
q from nearest-neighbors to essentially infinite, viz. values of q of (a) 2, (b) 3, (c) 7, (d) 15, (e) 25. 
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FIG. 4: (Color online) SMF and IMF models, with /3=2, for ranges of interaction q = 2, 15, 25. 



elastic and entropic contributions to the step-step inter- 
actions take the form v(x) = Ax~ 2 . In Refs. [40, 44 , it 
was shown that this kind of potential gives a good de- 
scription of the TWD for interacting steps. That is our 
motivation to extend our SMF model to this potential. 
We used the same approximation made in the previous 
section but now with this potential, finding 

p i0;q \s) =Ae~^ AC3H ^ s ~ i(3/2)A ^j=^ sHj ~ 1)r \ (51) 
The functional form of the SMF model is similar to 



the one of the zeroth-order model proposed in Ref. [40] . 
Fig. [6jb) shows that Eq. (51) accounts well for the nu- 
merical data. Remarkably, in spite of its simplicity, the 
SMF approximation describes pt '^ (s) very well even for 
large q. 



VI. CONCLUSION 

In many cases our knowledge about the statistical be- 
havior of a many-particle system is reduced to the spac- 
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ing distribution functions due to our inability to obtain 
numerically, experimentally, or analytically the joint dis- 
tribution function P/v(#i, ■ • * > x n', P)- Since the informa- 
tion contained in p^ q \s) is very limited, p^ n ' q \s) with 
n > should also be calculated to obtain additional phys- 
ical information about the system. 

In the IIA, the statistical behavior of the system is 
completely determined by p( 0;q \s). This is a crude ap- 
proximation which becomes exact only in the case of 
nearest-neighbor interactions. However, we showed that 
the IIA can be applied satisfactorily in many cases. One 
of these cases is Dyson's Brownian model for low values 
of q. In this case we expect that the correlations between 
gap sizes are weak, justifying the use of the IIA. Nat- 
urally, for large values of q the IIA gives poor results. 
We found also that the IIA — with the suitable choice of 
p( 0]q \s) — is a reasonable approximation for p( n;q \s) in 
the point-island model for epitaxial growth, at least in 
the ID and 2D. 

The EWS can be applied in systems where p( n;q \s) 
can be characterized by its behavior in the limits s«l 




and s ^> 1. This is the case of Dyson's Brownian model 
with q = 1 and q — N — 1 — >• oo. However the EWS does 
not give good results for this model in the case of 1 < q < 
N — 1. The EWS can also be used satisfactorily in the 
case of the 2D epitaxial growth mentioned previously. 

The mean-field approximation proposed in this paper 
is a general tool to study classical interacting particles 
with finite or full interaction range. The SMF and IMF 
models are more general approximations than the IIA or 
the EWS. In particular, the IMF model given by Eq. (|50|) 



can be interpreted as a generalization of the EWS for 
finite range of interaction. Additionally, the functional 
form of Eq. (44) can be used in different contexts [4T| |42], 
suggesting that this kind of model is related with other 
models for different systems. 

In order to generalize the SMF and IMF models for n > 
in Dyson's Brownian model, we calculated explicitly a n 
from the results given in Ref. pQ. We found that one of 
the effects of the finite range of interactions is to change 
a n . For n < q — 1, a n , has the same functional form 
of the case of complete range of interactions and it does 
not depend on q. For the case n > q — 1 we found a q 
dependence in a n . 

The IMF gives a good description of the Dyson's Brow- 
nian model at least for /3 = 1, 2 and 4. In spite of their 
simplicity, the SMF and IMF models give excellent re- 
sults for different kinds of interaction potentials. In par- 
ticular, we also tested the SMF approximation for the 
v(r) = A/s 2 potential of interacting steps. 
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This work was supported by the NSF-MRSEC at the 
University of Maryland, Grant No. DMR 05-20471 and 
DOE via University of Tennessee, with ancillary support 
from the Center for Nanophysics and Advanced Materials 
(CNAM). 



10 



Appendix A: Mean-field approximation 

As a first approximation we replace S n —> (S) = L/N 
for n ^ m. In this spirit, the average value over the 
ensemble of particles at time t of the quotients in Eq. ( [32] ) 
can be written as 



1 



— >> 



Ci 



and 



Co 



Sra + ECl 



— » 



S + (j-l)(S) : 



(Al) 



(A2) 



where Ci and C2 are constants. Naturally, these con- 
stants can be interpreted as a renormalization of the 
strength of the interaction between particles. 

We formulate an improved mean-field model by ana- 
lyzing more carefully the dependence on S m in Eq. (32). 
Following Refs. [22j [23], we write 



Appendix B: Exact behavior of the Dyson's 
Brownian model for s < 1 



From Eqs. (66) and (67) of Ref. pQ one can find the 
behavior of p( n;q \s) for small values of s. The case n < 
q — 1 is given by Eq. (66) of Ref. [1 . Since &i <C 

(S), we find 



POO POO 

P„+i(5i, • • • , S n+1 ) cc / dS n+2 ■■■ dS^aSMS), 
Jo Jo 

(Bl) 

where 

q-1 q-2 

as) = n s i n^+^+i)^ • • • - ( B2 ) 

i=l i=l 



and 



(B3) 



m+1 



■i) 



(A3) 



ft (5), (A4) 



where = Ci + 1) (S 2 )). Note that, in gen- 

eral, (S 2 ) is a function of t. However, in our mean-field 
approach the average of S 2 at time t is replaced by the av- 
erage in the stationary state (• • • ) st . In order to simplify 
our notation, in the text we omit the subscript indicating 
stationary state. In the same way we write 



1 



with q(S) = Ci/(j (S)). Finally, we use 



S ' 



Again, C\ and C2 are constants. 



ft (S), (A5) 



(A6) 



After integration, Eq. (Bl) reduces to 

n+l n 

P„+i(Si,"- ,S n+1 ) oc ]]s?l[(s i + s i+1 )t 



i=l i=l 

x(5i + --- 



■S n+ i)^(5)(B4) 



Finally, using Eq. (B4) with Eq. (12) we find 



/3(n + l)(rt + 2) 



p( n;q \s) oc s n +^?=o(^+ 1 - i ) oc g —'T'-'-' +n^ (B5) 
For n > g - 1 we use Eq. (67) of Ref. pQ. Then 

g— 1 n+l— j 

n n (B6) 

3 = 1 i=l 



Substituting Eq. (B6) into Eq. (12) leads to 



(^;<?) 



( s ) oc s /3El=o(n-i+l)+n = g § (3+2n-g)0+n (B?) 
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